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Abstract 



The Relativistic Random Phase Approximation (RRPA) is derived from the 
Time-dependent Relativistic Mean Field (TD RMF) theory in the limit of 
small amplitude oscillations. In the no-sea approximation of the RMF the- 
ory, the RRPA configuration space includes not only the ususal particle-hole 
p/i-states, but also a/i-configurations, i.e. pairs formed from occupied states 
in the Fermi sea and empty negative-energy states in the Dirac sea. The 
contribution of the negative energy states to the RRPA matrices is exam- 
ined in a schematic model, and the large effect of Dirac sea states on isoscalar 
strength distributions is illustrated for the giant monopole resonance in ^^^Sn. 
It is shown that, because the matrix elements of the time-like component 
of the vector meson fields which couple the a/i-configurations with the ph- 
configurations are strongly reduced with respect to the corresponding matrix 
elements of the isoscalar scalar meson field, the inclusion of states with unper- 
turbed energies more than 1.2 GeV below the Fermi energy has a pronounced 
effect on giant resonances with excitation energies in the MeV region. The 
influence of nuclear magnetism, i.e. the effect of the spatial components of 
the vector fields is examined, and the difference between the non-relativistic 
and relativistic RPA predictions for the nuclear matter compression modulus 
is explained. 
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I. INTRODUCTION 



Models based on the Relativistic Mean Field (RMF) approximation provide a microscopic 
self-consistent description of nuclear structure phenomena (for reviews, see Refs. In 
this framework classical equations of motion are derived self-consistently from a fully rela- 
tivistic Lagrangian. Vacuum polarization effects, as well as Fock exchange terms, are usually 
not taken into account explicitly. This framework, however, is based on an effective theory 
and the parameters of effective interactions are determined from a set of experimental data. 
In adjusting the parameters of the effective Lagrangian, a large part of vacuum polarization 
effects and effects of exchange terms are already taken into account. In fact, both contribu- 
tions have been treated explicitly in nuclear matter and in some finite spherical nuclei 
and it has been shown that these contributions are not small. However, by renormalizing the 
parameters of the effective Lagrangian, virtually identical results for ground-state properties 
have been obtained without the inclusion of vacuum polarization and exchange terms. Es- 
sential for a quantitative description of properties of complex nuclei are the non-linear terms 
in the meson sector |T^, which in a simple way include an effective density dependence of 
the meson coupling parameters. 

The relativistic mean-field models have been mostly applied in the description of ground- 
state properties of nuclei all over the periodic table. In several cases this framework has also 
been very successfully applied to excited states. For example, the cranked relativistic mean- 
field model |TT| describes a large variety of phenomena in rotational bands of superdeformed 



nuclei |jT2|-|14|. Another example is the Time-dependent RMF model which has been used 



to describe the dynamics of giant resonances in nuclei ITsHlTj. From the time evolution of 
multipole moments of the single-particle density, the excitation energies of giant resonances 
can be determined. Excellent agreement with experimental values of isoscalar and isovector 
giant monopole, giant quadrupole and the isovector giant dipole resonances has been ob- 
tained. The disadvantage of the time-dependent approach is that in some cases the large 
computational effort prevents an accurate description of excited states as, for instance, the 
dynamics of low-lying collective excitations. 

The Relativistic Random Phase Approximation (RRPA) represents the small amplitude 
limit of the time-dependent relativistic mean-field theory. Some of the earliest applications 
of the RRPA |]TB|-P^ to finite nuclei include the description of low-lying negative parity 
excitations in ^^O ||18[, and studies of isoscalar giant resonances in light and medium nu- 



clei |T^. These RRPA calculations, however, were based on the most simple, linear a — to 



relativistic mean field model. Only recently non-linear meson self-interaction terms have 
been taken included in RRPA calculations []23| , |2^ . However, in these calculations the RRPA 
configuration space included only ordinary particle-hole pairs. This seems a reasonable ap- 
proximation, since the states formed from occupied positive-energy states in the Fermi sea 
and empty negative energy states in the Dirac sea, have unperturbed energies more than 
1.2 GeV below the Fermi level. It turned out, however, that excitation energies of isoscalar 
resonances calculated in this way were very different from those obtained in the TD RMF 
approach with the same effective interactions p5|j26[| . 

Since it is well known that in the non-relativistic framework the RPA corresponds to the 
small amplitude limit of time-dependent Hartree-Fock (TDHF) (see, e.g., Ref. [0), these 
discrepancies remained an open puzzle for a couple of years. Only recently it has been shown 
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28| , |29[] that, in order to reproduce results of time-dependent relativistic mean-field calcu- 
lations for giant resonances, the RRPA configuration space must contain negative-energy 
states from the Dirac sea. In principle, if the Dirac sea were fully occupied, these config- 
urations would be forbidden by the Pauli principle. However, in the no-sea approximation 
the negative-energy states do not contribute to the nucleon densities, i.e. these states are 
not occupied. It is, thus, possible to form ah pairs (a empty state in the Dirac sea, h occu- 
pied state in the Fermi sea) and include them in the RRPA configuration space. Although 
formally possible, and also necessary in order to preserve symmetries, this procedure raises 
some serious conceptual problems because the ah configurations have negative unperturbed 
excitations energies. This means that the energy surface is no longer positive definite. The 
static solutions of the RMF equations correspond to saddle points on the energy surface, 
rather than to minima as in non-relativistic Hartree-Fock calculations. It is not, therefore, 
a priori clear that small amplitude oscillations around these stationary solutions will be 
stable. Furthermore, it is not obvious why configurations with unperturbed energies more 
than 1.2 GeV below the Fermi level have such a pronounced effect on the excitation energies 
of giant resonances in the MeV region. In non-relativistic calculations, for instance, it has 
been noted [0 that p/i-configurations with very large excitation energies affect the position 
of the spurious mode, but they have no effect on the excitation energies of giant resonances. 

The purpose of the present investigation is to clarify some of these problems and to 
understand in a better way the relation between the RRPA and the TD RMF in the no- 
sea approximation. The paper is organized as follows: the time-dependent RMF model 
is analyzed in Sec. II. In Sec. Ill the RRPA is derived from the TDRMF equations in 
the limit of small amplitude motion. In particular, we discuss the importance of a/i-pairs 
in the RRPA configuration space. In Sec. IV we introduce a relativistic extension of the 
Brown-Bolsterli model and solve it in the linear response approximation. It is shown how 
large matrix elements, which couple the a/i-sector with the p/i-configurations, can arise in 
the RRPA matrices. The large effect of Dirac sea states on isoscalar strength distributions 
is illustrated in Sec. V. The results are summarized in Sec. VI. 



II. THE TIME-DEPENDENT RELATIVISTIC MEAN-FIELD MODEL 

In quantum hadrodynamics the nucleus is described as a system of Dirac nucleons which 
interact through the exchange of virtual mesons and photons. The model is based on the 
one-boson exchange description of the nucleon-nucleon interaction. The Lagrangian density 
reads [|l| 

C = ij) {i'-) ■ d — ra) ip + -da ■ da — -ni^a'^ 

- i'[9a(y + 9l,1-^ + 9pl ■ pr + 67 ■ A ^— • (1) 

Vectors in isospin space are denoted by arrows, and bold-faced symbols will indicate vectors 
in ordinary three-dimensional space; a dot denotes the scalar product in Minkowski space 
{'y -uj = 'y'^ujfj^ = 7ot^o — ^^)- The Dirac spinor ip denotes the nucleon with mass m. rricr, m^, 
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and TUp are the masses of the cr-meson, the a;-meson, and the p-meson, and g^, guj, and gp are 
the corresponding couphng constants for the mesons to the nucleon, and /Tic = 1/137.036. 
Eq.(|l]), Q'^'^ , R^'^ , and F'^'^ denote the field tensors of the vector fields uj, p, and of the 
photon, respectively: 



(2) 



If the bare masses m, m^, and irip are used for the nucleons and the uj and p mesons, there 
are only four free model parameters: g„, g^ and gp. Their values can be adjusted to 
experimental data of just few spherical nuclei. This simple model, however, is not flexible 
enough for a quantitative description of properties of complex nuclear systems. An effective 
density dependence has been introduced by replacing the quadratic cr-potential \m?^a'^ 
with a quartic potential U (a) 



U{<y) 



1 2 2 Q 1 A 

+ ^92CT + -gsa 



(3) 



The potential includes the nonlinear a self-interaction, with two additional parameters g2 
and (?3. The corresponding Klein-Gordon equation becomes nonlinear, with a a-dependent 
mass m'^{a) = + g2(T + g^a^. More details on the relativistic mean- field formalism can 
be found in Refs. [|l]-|] . 

From the Lagrangian density the set of coupled equations of motion is derived. The 
Dirac equation for the nucleons 



{aHV - V(r, t)]+V{v, t) + /?(m - ^(r, t))} ^i. 



(4) 



If one neglects retardation effects for the meson fields, the time- dependent mean-field po- 
tentials 



S{r,t) = -g^a{r,t) , 

Vpir, t) = g^ujpiv, t) + gpTppijc, t) + eA^(r, t) 



(5) 



are calculated at each step in time from the solution of the stationary Klein-Gordon equations 
-A + m^] a{T,t) = -g^ Ps(r,t) - 5-2 cr^(r,t) - 93 cr^(r,it) , 

Pfi{r,t) = gp j^{r,t) , 



-A - 
-A 



-AAp{r,t) = ejcp{r,t) . 



(6) 



This approximation is justified by the large masses in the meson propagators. Retardation 
effects can be neglected because of the short range of the corresponding meson exchange 
forces. In the mean-field approximation only the motion of the nucleons is quantized, the 
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meson degrees of freedom are described by classical fields which are defined by the nucleon 
densities and currents. The single-particle spinors {i = 1,2,..., A) form the A-particle 
Slater determinant The nucleons move independently in the classical meson fields, 

i.e. residual two-body correlations are not included, and the many-nucleon wave function is 
a Slater determinant at all times. The sources of the fields in the Klein- Gordon equations 
are the nucleon densities and currents calculated in the no-sea approximation 

A 

i=l 
A 

i=l 
A _ 

i=l 
Z 

j,^{r,t)^^^,{r,t)j^^,{r,t) . (7) 

i=l 

where the summation is over all occupied states in the Fermi sea. In the no-sea approxima- 
tion the negative-energy states do not contribute to the densities and currents, i.e. vacuum 
polarization is explicitly neglected. However, as already discussed in the Introduction, this 
is an effective theory with the parameters of the Lagrangian determined from a set of experi- 
mental data. In adjusting the parameters of the effective Lagrangian, a large part of vacuum 
polarization effects is therefore already taken into account. It should be emphasized that 
the no-sea approximation is essential for practical applications of the relativistic mean-field 
model. 

The stationary solutions of the relativistic mean-field equations describe the ground- 
state of a nucleus. They correspond to stationary points on the relativistic energy surface. 
The Dirac sea, i.e. the negative energy eigenvectors of the Dirac hamiltonian, is different 
for different nuclei. This means that it depends on the specific solution of the set of non- 
linear RMF equations. The Dirac spinors which describe the ground-state of a finite nucleus 
(positive energy states) can be expanded, for instance, in terms of vacuum solutions, which 
form a complete set of plane wave functions in spinor space. This set is only complete, 
however, if in addition to the positive energy states, it also contains the states with negative 
energies, in this case the Dirac sea of the vacuum. Positive energy solutions of the RMF 
equations in a finite nucleus automatically contain vacuum components with negative energy. 
In the same way, solutions which describe excited states, as for instance states with different 
angular momenta which are solutions of the cranked RMF equations, contain negative energy 
components which correspond to the ground-state solution. 

This is also true for the solutions of the time-dependent problem. Although for the 
stationary solutions the negative-energy states do not contribute to the densities in the no-sea 
approximation, their contribution is implicitly included in the time-dependent calculation. 
The coupled system of RMF equations describes the time-evolution of A nucleons in the 
effective mean-field potential. Starting from the self-consistent solution which describes the 
ground-state of the nucleus, initial conditions can be defined which correspond, for instance, 
to excitations of giant resonances in experiments with electromagnetic or hadron probes. 
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For example, the one-body proton and neutron densities can be initially deformed and/or 
given some initial velocities. The resulting mean-field dynamics can be described by the 
time-evolution of the collective variables. In coordinate space for example, these will be the 
multipole moments of the density distributions. At each time t, the Dirac spinors ijji{t) can 
be expanded in terms of the complete set of solutions of the stationary Dirac equation 

^.(r,t) = 5:cfc(t)4°V)e-^^'^-* + Eca(t)^i°Hr)e-^-* , (8) 

k a 

where the index k runs over all positive energy eigen-solutions Ek > (hole states h in the 
Fermi sea, and particle states p above the Fermi sea), and the index a denotes eigen-solutions 
with negative energy < 0. We follow the time evolution of A Dirac spinors which at time 
t = form the Fermi sea of the stationary solution. This means that at each time we 
have a local Fermi sea of A time-dependent spinors which, of course, contain components of 
negative-energy solutions of the stationary Dirac equation. One could also start with the 
infinitely many negative energy solutions ipai^^t = 0) {Sa < 0), and propagate them in time 



with the same hamiltonian h(t). Since the time-evolution operator is unitary |27 



tdt {^^\^a) = {^i\h^ - h\^a) = 0, (9) 

the states which form the local Dirac sea are orthogonal to the local Fermi sea at each time. 
This is the meaning of the no-sea approximation in the time-dependent problem. For small- 
amplitude oscillations around the stationary solution, the coefficients Ca{t) of the negative 
energy components in (H) are, of course, small. 

We will first consider linear relativistic mean-field models ((72 = fi's = in (^). In the 
instantaneous approximation, i.e. neglecting the time derivatives in the Klein-Gordon 
equations, the solutions for the mean-fields are calculated from 

cr{r,t) = 9a J D^{r,r')ps{r',t)d^r , 

w^(r,t) = J D^{r,v')j^{r',t)dh- , 

= gp J Dp{r,r%{r',t)d^r , 

A^(r, t)=ej Dphotonir, r')j,p{r', t)d^r . (10) 

The propagators have the Yukawa form 

1 p-m</.|r-r'| 

where (p denotes the mesons a, u, p, and the photon. The plus (minus) sign is for vector 
(scalar) fields. In the non-linear case an analytic solution of the Klein-Gordon equation is, 
of course, no longer possible. The corresponding meson field is a non-linear functional of 
the density and currents. 

The relativistic single-particle density matrix reads 

p(r,r',t)=X:iV'.(r,t))(^,(r',t)| . (12) 

i=l 
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If the Dirac spinor is written in terms of large and small components 

l^.(r,t)) = 



igi{r,t) 



(13) 



the density matrix takes the form 



p(r,r',t) 



E/.(r,t)/;(r',t) -zE/.(r,tW(r',t)^ 

i=l i=l 
i=l i=l 



(14) 



Further, a relativistic two-body interaction is defined 

V = [ rfVirf\2V^^(ri)V^^(r2)F(ri,r2)^(ri)^(r2) , 



(15) 

where '?/'^and ip are the Dirac field creation and annihilation operators, and 

r(ri, r2) = D.(ri, r^) P'-'^P^'^ + D^{t^, v^) (l - cx^'^cx^^A . (16) 



In order to simplify the notation, we omit the p-meson and the photon, though they are, of 
course, included in actual applications of the relativistic mean-field model. Their contribu- 
tion to the matrix elements of V^(ri, r2) is, however, much smaller than that of the a and uj 
mesons. 

By introducing an arbitrary complete spinor basis (the indices k,l, . . . denote both pos- 
itive and negative energy states), the two-body interaction operator can be written in the 
form 



^ kk'W 



(17) 



The single-particle equation of motion corresponds to the time-dependent relativistic Hartree 
problem 



with the Dirac Hamiltonian 



and the mass operator 



h{p) = CK p + (3{m + S(p)), 



^ki{p) — '^VkVik'Pk'i' ■ 

k'V 



The corresponding equation of motion for the density operator reads 

idtp = [h{p)^p\ ' 

in full analogy with the non-relativistic Hartree- Fock problem (see, e.g., Ref. 



(18) 



(19) 



(20) 



(21) 



7 



In expressing the TD RMF equations (^y) in terms of a relativistic two-body interaction, 
we have ehminated the meson degrees of freedom by using the Yukawa form (p!lD of the 
meson propagators. This apphes, of course, only to Lagrangians that do not contain non- 
hnear meson self-interactions. The non-linear couplings are, however, essential for a realistic 
description of nuclear properties. Formally, also in this case the Klein-Gordon equations 
can be solved at each step in time, and the resulting meson fields are non-linear functional 
of the densities and currents. The Dirac operator has still the form of Eq. (|T^), but the 
mass-operator ^^^(p) becomes a much more complicated functional of the single-particle 
density. 

The numerical solution of the full time- dependent problem with non-linear meson self- 
interactions does not present particular difficulties (see Refs. [ITsH TTH). Much more difficult, 
however, is to eliminate the meson degrees of freedom and to derive a relativistic two-body 
interaction in the general case of large amplitude motion. This has only been done in 
the small amplitude limit [^]. The cx-field and the scalar density ps are expanded in the 
neighborhood of the stationary ground-state solutions 



a(r,t) = (T^"^(r)+(5(T(r,t) 



(22) 
(23) 



The corresponding Klein- Gordon equation 
up to terms linear in 5a we obtain 



for the (T-field is solved by linearization, i.e. 



-A + m^(r) 5 cr(r,t) = -g„5 Psir,t) 



(24) 



with 



■ml(r) 



(T=a{0)(r) 

The cr-meson propagator is defined by the equation 



(25) 



-A + m^(r) 



DJr,r') 



-(5(r 



(26) 



23 



In 



and it has been determined numerically in the RRPA calculations of Refs. 
the following only the small amplitude limit will be studied, and therefore we do not need 
to worry about the more general problem of large amplitude motion. 



III. THE SMALL AMPLITUDE LIMIT OF TD RMF AND THE RELATIVISTIC 

RPA 

In this section we study the response of the density matrix p{t) to an external one-body 
field 

F{t) = Fe"*"* + h.c. , (27) 

which oscillates with a small amplitude. Assuming that in the single-particle space this field 
is represented by the operator 
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kl 



the equation of motion for the density operator is 

idtp = [h{p) + f{t),p\ , 
In the hnear approximation the density matrix is expanded 

m = p(°) + 5m , 



(28) 



(29) 



(30) 



where p^'^^ is the stationary ground-state density. From the definition of the density matrix 
(|1^), it follows that p(t) is a projector at all times, i.e. p(t)^ = p(t). In particular, this 
means that the eigenvalues of p^^^ are and 1. In the non-relativistic case particle states 
above the Fermi level correspond to the eigenvalue 0, and hole states in the Fermi sea 
correspond to the eigenvalue 1. In the relativistic case, one also has to take into account 
states from the Dirac sea. In the no-sea approximation these states are not occupied, i.e. 
they correspond to the eigenvalue of the density matrix. We will work in the basis which 
diagonalizes 



Pkl - 5klPk 



(0) 



for unoccupied states above the Fermi level (index p) 

1 for occupied states in the Fermi sea (index h) 

for unoccupied states in the Dirac sea (index a) 



Since p{t) is a projector at all times, in linear order 

p^^'^5p + 5pp^^'^ = 6p 



(31) 



(32) 



This means that the non-vanishing matrix elements of 6p are: Spph, Sphp, Spah, and Spha- 
These are determined by the solution of the TD RMF equation (P5|). In the linear approxi- 
mation the equation of motion reduces to 



idt5p= \U^\Sp 



+ 



9h .(0) 
d-p^P^P 



+ 



(0) 



where 



dh dh dh 

9p ^ 9pph dphp 



„ dh dh 

1^ ~ "Pah + — Opha 



ah ^P»h 



dp, 



'ha 



(33) 



(34) 



In the small amplitude limit 5p will, of course, also display a harmonic time dependence 
g-iaji_ 'jaking into account the fact that hfi = 6kitk is diagonal in the stationary basis, we 
obtain 



(a; — ep + eh)Spph — fph + ^ Vph'hp'^Pp'h' + Vpp'hh'^Ph'p' + X! ^ph'ha'^Pa'h' + Vpa'hh'^Ph'a' 

p'h' a'h' 

fah + X! ^ah'hp'SPp'h' + Vap'hh'^ Ph'p' + X! ^ah'ha'^ Pa'h' + Vaa' hh'^ Ph'a' 
p'h' a'h' 

fhp + X! ^hh'pp'SPp'h' + Vhp'ph'SPh'p' + X! ^hh'pa'SPa'h' + Vha'ph'^ Ph'a' 
p'h' a'h' 

{lO — eh + ^a)^Pha = fha + X! ^hh'ap'^ Pp'h' + Vhp' ah'^ Ph'p' + ^hh' aa' ^ Pa' h' + Vha' ah'^ Ph' a' (35) 

p'h' a'h' 



(CJ - e„ + ehjSpah 

- eh + ep)6php 
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or, in matrix form 



1 



-1 J 

The RRPA matrices A and B read 



A B 

B* A* 



X 
Y 



F 
F 



A 
B 



^h)^aa'^hh' 



+ 



ph'hp' ^ph'ha' 



Vr 



ah'hp' Vah'ha' 



Vpp'hh' Vpa'hh' 
Vap'th' Vaa'hh' 



(36) 

(37) 
(38) 



and the amphtudes X and Y are defined 



X = ( ^Ppf' \ Y = ( ^^^^ \ 



(39) 



The vectors which represent the external field contain the matrix elements 

F = {p'^ V F - 

J ah 



fhp 
fha 



(40) 



In conventional linear response theory (see, e.g., Ref. [^) the polarization function Ilpqp'q'{uj) 
is defined by the response of the density matrix to an external field with a harmonic time 
dependence 



^Ppq — '^^pqp'q'i^) fp'q' ■ 
p'q' 

Its spectral representation reads 

(oiv^U i/i)(/i|4v^,|0) (o|4v^.|/i)(/i|v^t^ 10) 



(41) 



00 - + Eo + ir] 



u + Ef, - Eo + ir] 



(42) 



where the index n runs over all excited states \n) with energy E^^. In the RPA approximation 
the polarization function is obtained by inverting the matrix 



Il(uj) 



uj + iT] \ ( A B 
-u -iri ) ~ \B* A* 



n -1 



n(a;) is the solution of the linearized Bethe-Salpeter equation 

n(cj) = n°(tu) + ii^{uj)vii{u) , 

where the free polarization function is given by 



io - tk + ei + trj 

The eigenmodes of the system are determined by the RPA equation 



(43) 



(44) 



(45) 
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(46) 

In principle, this is a non-Hermitian eigenvalue problem. In the non relativistic case, how- 
ever, it can be reduced to a Hermitian problem of half dimension, if the RPA matrices are 
real and if {A + B) is positive definite. In this case one can also show that the eigenvalues 
Q't are positive, i.e., the RPA eigenfrequencies are real (see |27||). 



The relativistic case is much more complicated. From Eq. (|37| ) we notice that the matrix 
{A + B) is not positive definite. The ah configurations have large negative diagonal matrix 
elements eah = — c/i < —1.2 GeV, and the RRPA equation can no longer be reduced 
to a Hermitian problem of half dimension. In this case it is also not clear whether the 
eigenfrequencies are necessarily real, because the stability matrix 



•5= " rJ (47) 




is no longer positive definite. Rather than minima, the solutions of the RMF equations are 
saddle points ^2|] in the multi-dimensional energy surface, and the Thouless theorem ^3 



which states that a positive definite stability matrix S leads to a stable RPA equation with 
real frequencies, does not apply. 

However, the opposite is not true: if the stability matrix is not positive definite, it does 
not automatically follow that the eigenvalues of the corresponding RPA matrix are not real. 
In fact, cases like this occur also in the non relativistic RPA in the neighborhood of phase 
transitions, where the interaction V is very large and attractive. The positive energies 
£p — £h on the diagonal of the stability matrix are not large enough, as compared to the 
matrix elements of V, to guarantee positive eigenvalues of S. In the relativistic case the 
energies on the diagonal Ea — Sh are negative. Even for small matrix elements of V the 
stability matrix S will have negative eigenvalues. However, as long as the diagonal part 
dominates, i.e. as long as we are not in a neighborhood of a phase transition, the RRPA 
eigenfrequencies are real. This can be easily demonstrated if instead of the RPA amplitudes 
X and Y , we define the generalized coordinates Q and momenta P 

g = i=(X-n, P=^(X + n- (48) 

In the small amplitude limit the time-dependent mean field equations take the form of 
classical Hamiltonian equations (for details see Ref. |2^, Chapt. 12) for the Hamiltonian 
function 

nP.Q) = \[P'^ -P)^''[lp'^)+\{Q' Q)s[q*)^ (49) 
with the inertia tensor 




M=[ % f:] . (50) 



The large negative diagonal matrix elements are also present in the inertia tensor. If the 
off-diagonal matrix elements are not too large, a negative inertia and a negative curvature 
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will again result in real frequencies. In all applications of RRPA we have found real frequen- 
cies, though in none of these cases the stability matrix S was positive definite. This also 
explain why the time- dependent RMF equations have stable solutions which describe oscil- 
lations with real frequencies around the static solution, although the static solution itself 
corresponds to a saddle point. 

The solution of the RPA equations in configuration space is much more complicated in 
the relativistic case. Firstly, because in addition to the usual p/i-states, the configuration 
space includes a large number of a/i-states. A further complication arises because the full 
non-Hermitian RPA matrix has to be diagonalized, even in cases when the matrix elements 



are real. The usual method ||2^ , which reduces the dimension of the RPA equations by half 
does not apply. 

Summarizing the results of this section, we have shown that the relativistic RPA repre- 
sents the small amplitude limit of the time- dependent RMF theory. However, because the 
RMF theory is based on the no-sea approximation, the RRPA configuration space includes 
not only the ususal ph-states, but also a/i-configurations, i.e. pairs formed from occupied 
states in the Fermi sea and empty negative-energy states in the Dirac sea. At each time t 
the occupied positive energy states can have non-vanishing overlap with both positive and 
negative energy solutions calculated at t = 0. If the density matrix p{t) is represented in the 
basis which diagonalizes the static solution p^^\ it contains not only the usual components 
6pph with a particle above the Fermi level and a hole in the Fermi sea, but also components 
Spah with a particle in the Dirac sea and a hole in the Fermi sea. 

One of the important advantages of using the time-dependent variational approach is 
that it conserves symmetries. It is well known from non-relativistic time-dependent mean 
field theory that symmetries are connected with zero energy solutions of the RPA, i.e. the 
Goldstone modes, and it is one of the advantages of RPA that it restores the symmetries 
broken by the mean field. This has already been realized in the early studies of symmetry 
conservation in RRPA, and it has been emphasized by Dawson and Furnstahl in Ref. pi 



that it is essential to include the ah configuration space in order to bring the Goldstone 
modes to zero energy. 

However, it was not anticipated that negative energy states in the RRPA configuration 
space could have dramatic effects on the excitation energies of giant resonances, as we will 
show in the following sections. It is not obvious that basis states with unperturbed energies 
more than 1.2 GeV below the Fermi energy, can have a big influence on giant resonances 
with excitation energies in the MeV region. In the following section we will study a simple 
model which provides a deeper insight into this problem. 



IV. A SEPARABLE MODEL 

The model studied in this section represents a relativistic extension of the Brown and 



Bolsterli model [0, which has played an essential role in the understanding of the micro- 
scopic picture of collective excitations. 

The single-particle basis consists of 4 states, each of them fi-fold degenerate (i^ = 1, . . . 
The first two states (1 and 2) correspond to particle levels with the free mass m, the states 3 
and 4 correspond to the negative energy levels with free mass — m. The model Hamiltonian 
reads 
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H = Ho- \,S^S + X.V^V , 
where Hq is the Hamiltonian which describes free Dirac particles 

A ^ 



(51) 



(52) 



with 
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(53) 



In Eq. (|5^) mj = m is the free mass of particle i, Pi = p denotes its momentum, and 
Ei = Eq <^ m induces a small splitting between the levels 1 and 2 (and, of course, between 3 
and 4). 

The interaction consists of an attractive scalar field S and a repulsive vector field V 



A 
i=l 



fO 1 \ 
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Vo -1 y 
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^0 1 ON^ 
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(54) 



with the strength parameters Xg and A^. In the formalism of second quantization the oper- 
ators -f^o, "S* and V take the forms 

V 

+m c\uCiu + etc. 



2v C^yC^v Ca„C, 



I ^0 v;^ t _ t it 

^ u 

S = ^ c\^C2y — cl^c^u + h.c. , 
^ = Y1 ^iy^2u + cLc4^ + h.c. 



(55) 
(56) 
(57) 



At the mean field level the diagonalization of the Dirac operator 



result in the eigenvalues 

.. = £ + |. 



/m + f 
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(58) 
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and the corresponding eigenvectors are 





9 

VO/ 



f 



\9J 



f-3\ 


/ 

V / 



respectively. We use the notation 

E = ^Jp^ + m^,, 

with 



/ = cos-, 



tan — 



P 



2 m + E 
A reahstic choice of single particle energies is 



( \ 

-9 


V / J 



. 

^ = sm - , 



eph = ep - = £0 ~ 10 MeV , 

eah = ea-eh = -2E + - -2 GeV , 

ea'h = ea' -eh = -2E ~ -2 GeV . 



(60) 



(61) 



(62) 



(63) 



In realistic calculations the ratio between large and small components of the Dirac spinors 
is approximately //(? ~ 30, i.e., ~ 33°. In the basis (^) the matrices of the operators 5" 
and V are 



( 

COS0 



V — sin ( 



COS</) 



— sin I 






- sine 


COS( 



sin \ 




cos< 




(64) 
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We notice the essential matrix elements 



Sph = cos (f) Vph = 
Sah = - sin Vah -- 

Sa'h = Va'h 



1 



= . 



(65) 



(66) 



In analogy to the Brown-Bolsterli model, the unperturbed polarization function is 



[00 



'^^oP'FphFph 



+ 



uj^ - ef^ + IT] uj^ - e^f^ + 17] 
where the operators F, F' G {5*, V}. In particular. 



(67) 
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Kvi-^) = , . , (69) 

UJ^ — Eq + IT] 

Ksi^) = Ksi^) = o ° , ■ (70) 
The RRPA frequencies are the roots of the determinant 

The essential difference with respect to the non-relativistic Brown-Bolsterh model is the 
additional term 2ii^nsin^ — E"^ + irf) in the scalar polarization n^_g(u;). Without this 

term (i.e. = 0), the eigenfrequencies would be determined by 

uj^ = el-~ 2e^n{K - K) , (72) 

with the usual cancellation of scalar and vector interactions. With the additional term, 
states with unperturbed energies at ~ —2E are included in the RPA configuration space. 
The interaction between these states and the p/i-states is determined by the matrix elements 
of the scalar interaction 

Vah'hp = -K COS sin . (73) 

These matrix elements are not reduced by a similar term coming from the vector interaction. 
In our simplified model this vector-induced term vanishes as a consequence of the relativistic 
structure of the equations: while for a state from the Fermi sea the large component is the 
upper component of the spinor, a state from the Dirac sea has a large lower component of 
the spinor. Due to the 7-matrix structure of the vertex, the matrix elements of the vector 
interaction vanish. In realistic calculations these matrix elements do not vanish identically, 
but they are reduced by an order of magnitude as compared to the corresponding scalar 
terms. 

At excitation energies in the MeV region, uj <^ E in the denominator of the second term 
of Eq. ([68|) , and we obtain an energy-independent term {2flXsSm'^ (j))/E in the dispersion 
relation. The eigenfrequencies are now determined by 

u;'' = el- e,2niXs - A,) + e,2nXs sin^ ^ f , , (74) 

E + 2\lXs sm 

i.e. we find an additional repulsion for the collective state. 



V. AN ILLUSTRATIVE CASE: THE GIANT MONOPOLE RESONANCE 

The RPA equations can be solved either by diagonalizing the RPA matrix in configura- 
tion space (see Eq. (|46|) ), or the response function can be calculated by solving the Bethe- 
Salpeter equation (|^) in momentum space ||T9| , p3| . In both cases, of course, one first has to 
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determine the single-nucleon spinors and the mean-fields which correspond to the stationary 
solution for the ground-state. The Dirac-Hartree equations and the equations for the meson 
fields are solved self-consistently in the mean-field approximation. The eigenvalue problem 
is solved, for instance, by diagonalization in a spherically symmetric harmonic oscillator 
basis From the spectrum of single-nucleon states the RPA configuration space is built: 
particle-hole (ph) and antiparticle-hole (ah) pairs which obey the selection rules for angu- 
lar momentum, parity and isospin. The number of basis states is also determined by two 
cut-off parameters: the maximal ph-energy (e.^ — < -Emax) and the minimal ah-eneigy 
i^a — ^i> -S'min)- With this basis the RPA matrix is calculated for the same effective interac- 
tion that determines the ground-state, or the free polarization function 11'^'') is calculated in 
the response function method. Both methods require that the single-particle continuum is 
discretized. In order to smooth out the RPA strength function, the discrete strength distri- 
bution is folded by a Lorentzian of width F. In the response function method the folding is 
automatic if a finite value parameter iT is used in the denominators of Eqs.(^^-^), instead 
of the infinitesimal parameter if]. We have verified that identical results are obtained with 
both methods. 

The large effect of Dirac sea states on isoscalar strength distributions is illustrated in 
Fig. 1, where we display the isoscalar monopole RRPA strength in ^^^Sn calculated with 



the NL3 effective interaction [361 and the width of the Lorentzian is F = 2 MeV. Recent 



experimental data are available for the isoscalar giant monopole resonance in ^^^Sn [|35 



The solid curve represents the full RRPA strength and it displays a pronounced peak at 16 
MeV, in excellent agreement with the measured value of 15.9 MeV Giant monopole 

resonances in spherical nuclei are in best agreement with experimental data, when calculated 
with effective Lagrangians with a nuclear matter compression modulus in the range 250-270 
MeV |0123J23]. The nuclear matter incompressibility of the NL3 effective interaction is 272 



MeV. 

The long-dashed curve in Fig. 1 corresponds to the to the case with no ah pairs in 
the RRPA configuration space. We notice that, without the contribution from Dirac sea 
states, the strength distribution is shifted to lower energy. The position of the peak is 
shifted from ^ 16 MeV to below 10 MeV if ah pairs are not included in the RRPA basis. 
Quantitatively similar results are also obtained with other effective interactions. In Fig. 1 
we have also separated the contributions of vector and scalar mesons to the ah matrix 
elements. The dash-dot-dot (dash-dot) curve corresponds to calculations in which only 
vector mesons (scalar mesons) were included in the coupling between the Fermi sea and 
Dirac sea states. Both interactions were included in the positive energy particle-hole matrix 
elements. The resulting strength distributions nicely illustrate the dominant contribution 
of the isoscalar scalar sigma meson to the ah matrix elements, in complete agreement with 
the result obtained in the previous section for the schematic Brown-Bolsterli model. 

It is also interesting to examine the effect of the coupling via the spatial components of 
the vector meson fields, i.e. the term —a^^^a.^'^^ in the interaction of Eq. ([T6|) . In time- 
dependent calculations this coupling results from the nucleon currents. In Fig. 2 we display 
the isoscalar monopole RRPA strength in ^^^Sn calculated as follows: a) full RRPA (solid 
curve); b) without the matrix elements of the spatial components of the vector meson fields 
(dot-dashed curve); c) without the contribution of the Dirac sea to the matrix elements of 
the spatial components of the vector meson fields (dashed curve); and d) the free Hartree 
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response function. The currents do not contribute to the static polarizabihty and to the M_i 
moment. At finite frequencies, however, their contribution is attractive and it lowers the 
ISGMR energy by ~ 2 MeV.. Therefore, if the contribution of the spatial components of the 
vector fields is neglected, a better agreement with experimental values would be obtained 
with a lower nuclear matter incompressibility: K^o — 230 MeV. Incidentally, this lower 
value for the nuclear matter incompressibility is the one advocated by non-relativistic RPA 
calculations [^^^. It has already been noted in time-dependent RMF calculations []T^, as 



well as in recent relativistic RPA studies |29], that effective interactions which reproduce 



the IS GMR excitation energies in finite nuclei have a somewhat higher nuclear matter in- 
compressibility than the corresponding non-relativistic Skyrme or Gogny interactions. Here 
we point to a possible solution to this puzzle: the current terms in the matrix elements of 
the particle-hole interaction (^,^) are given by 




{pmkr)[c.Y,{v)]j=,\\h) = {plMkr) { " ^^'=' ] \\h) , (75) 



which is a typical relativistic term because it couples large and small components of a Dirac 
spinor. Since they change parity, terms of the type [(jYi]j=o cannot contribute to the giant 
monopole resonance in a non-relativistic calculation. 



VI. CONCLUSIONS 

In the last couple of years, several discrepancies have been reported between the results 
obtained with the Relativistic Random Phase approximation and the Time-Dependent Rela- 
tivistic Mean Field theory, when applied to the the description of small amplitude collective 
motion in atomic nuclei. 

In order to resolve this puzzle, in the present work we have derived the RRPA from 
the TDRMF equations in the limit of small amplitude motion. The relativistic single par- 
ticle density matrix p{t) has been expanded in terms of the stationary solutions of the 
ground-state. We have shown that the no-sea approximation, which is essential for practical 
application of the RMF theory in finite nuclei, leads to a fundamental difference between 
the relativistic and non-relativistic approaches. While in the non-relativistic case the time- 
dependent variation of the density 5p{t) = p{t) — p^^^ has only p/i-matrix elements (particle 
(p) above the Fermi surface, hole {h) in the Fermi sea), in the relativistic case 6p contains 
also a/i-matrix elements, where a denotes unoccupied states in the Dirac sea. The fact 
that states in the Dirac sea can be occupied is a direct consequence of the no-sea approxi- 
mation. In constructing the matrix Sp one has to take into account that a complete basis 
of single particle states contains both positive and negative energy solutions of the Dirac 



equation. Already in Ref. it has been shown that an RRPA calculation, consistent with 
the mean-field model in the no — sea approximation, necessitates configuration spaces that 
include both particle-hole pairs and pairs formed from occupied states and negative-energy 
states. The contributions from configurations built from occupied positive-energy states 
and negative-energy states are essential for current conservation and the decoupling of the 
spurious state. 

What is less obvious, however, is that the inclusion of negative-energy single particle 
states in the RRPA configuration space has such a dramatic effect on the calculated excita- 
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tion energies of isoscalar giant resonances. In a schematic model we have shown that, due to 
the relativistic structure of the RPA equations, the matrix elements of the time-hke compo- 
nent of the vector meson fields which couple the a/i-configurations with the p/i-configurations 
vanish. In realistic calculations these matrix elements do not vanish identically, but they are 
strongly reduced with respect to the corresponding matrix elements of the isoscalar scalar 
meson field. As a result, the well known cancellation between the contributions of the a and 
10 fields, which, for instance, leads to ground-state solution, docs not take place and we find 
large matrix elements coupling the a/i-sector with the p/i-configurations. In addition, the 
number of a/i-configurations which can couple to the p/i-configurations in the neighborhood 
of the Fermi surface is much larger than the number of p/i-configurations. This can increase 
the effect by enhancing the collectivity of the a/i-configuration space. 

The large effect of Dirac sea states on isoscalar strength distributions has been illustrated 
for the giant monopole resonance in ^^^Sn. We have also shown that currents cannot be ne- 
glected in the calculation of giant resonances. Of course they do not occur in the static case, 
i.e. in the calculations of the static polarizability or the M_i moment. At finite frequencies, 
however, time reversal invariance is broken and spatial components of the vector meson fields 
play an important role. This effect is known as nuclear magnetism. It is a genuine relativis- 
tic effect, because the matrix elements couple the large and small components of a Dirac 
spinor. Since the spatial components of the vector fields have the form —oS^^ol''^^ D^^{ri, r2), 
they result in an attractive contribution which lowers the value of the calculated excitation 
energies of giant resonances. This explains the difference between the non-relativistic and 
relativistic RPA results for the isoscalar giant monopole resonances in spherical nuclei, and 
the corresponding predictions for the nuclear matter compression modulus. 

ACKNOWLEDGMENTS 

P.R. acknowledges the support and the hospitality extended to him during his stay at the 
IPN-Orsay, where a large part of this work was completed. This work has been supported in 
part by the Bundesministerium fiir Bildung und Forschung under the project 06 TM 979 and 
by the Deutsche Forschungsgemeinschaft. It was also partially supported by the National 
Natural Science Foundation of China under grant No. 19847002, 19835010-10075080 and 
Major State Basic Research Development Program under contract No. G200077407. 



18 



REFERENCES 



B.D. Serot and J.D. Walecka, Adv. Nucl. Phys. 16, 1 (1986). 
P.-G. Reinhard, Rep. Progr. Phys. 52, 439 (1989). 

B. D. Serot, Rep. Prog. Phys. 55, 1855 (1992). 

P. Ring, Progr. Part. Nucl. Phys. 37, 193 (1996). 

C. J. Horowitz and B.D. Serot, Nucl. Phys. A399, 529 (1983). 

A. Bouyssy, J.F. Mathiot, N. Van Giai, and S. Marcos, Phys. Rev. C36, 380 (1987). 
P. Bernados, V.N. Fomenko, N. Van Giai, M.L. Quelle, S. Marcos, R. Niembro, and 
L.N. Savushkin, Phys. Rev. C48, 2665 (1993). 

C. J. Horowitz and B.D. Serot, Phys. Lett. 140B, 181 (1984). 

D. A. Wasson, Nucl. Phys. A535, 456 (1991). 

J. Boguta and A. R. Bodmer, Nucl. Phys. A292, 413 (1977). 

W. Koepf and P. Ring, Nucl. Phys. A493, 61 (1989). 

A.V. Afanasjev, J. K5nig, and P. Ring, Nucl. Phys. A608, 107 (1996). 

A.V. Afanasjev, G.A. Lalazissis, and P. Ring, Nucl. Phys. A634, 395 (1998). 

A. V. Afanasjev, J. Konig, and P. Ring, Phys. Rev. C60, 051303 (1999). 
D. Vretenar, H. Berghammer, and P. Ring, Nucl. Phys. A581, 679 (1995). 

B. Podobnik, D. Vretenar, and P. Ring, Z. Phys. A354, 375 (1996). 

D. Vretenar, G.A. Lalazissis, R. Behnsch, W. Poschl, and P. Ring, Nucl. Phys. A621, 
853 (1997). 

R.J. Furnstahl, Phys. Lett. 152B, 313 (1985). 

M. L'Huillier and N. Van Giai, Phys. Rev. C39, 2022 (1989). 

P.G. Blunden and M. McCorquodale, Phys. Rev. CSS, 1861 (1988). 

J.F. Dawson and R.J. Furnstahl, Phys. Rev. C42, 2009 (1990). 

C. J. Horowitz and J. Piekarewicz, Nucl. Phys. A511, 461 (1990). 

Z.Y. Ma, N. Van Giai, H. Toki, and M. L'Huillier, Phys. Rev. C42, 2385 (1997). 

Z.Y. Ma, H. Toki, and N. Van Giai, Nucl. Phys. A627, 1 (1997). 

N. Van Giai, Z.Y. Ma, H. Toki, and B.Q. Chen, Nucl. Phys. A649, 37c (1999). 

D. Vretenar, P. Ring, G.A. Lalazissis, and N. Paar, Nucl. Phys. A649, 29c (1999). 

P. Ring and P. Schuck, The Nuclear Many-Body Problem, Springer Verlag, New York 
1980. 

D. Vretenar, A. Wandelt, and P. Ring, Phys. Lett. B487, 334 (2000). 

Z.Y. Ma, N. Van Giai, A. Wandelt, D. Vretenar, and P. Ring, Nucl. Phys. A 685 (2001) 

(in press). 

P. Ring and J. Spcth, Nucl. Phys. A235, 315 (1974). 

G.E. Brown and M. Bolstcrh, Phys. Rev. Lett. 3, 472 (1959). 

T. Kohmura, Y. Miyama, T. Nagai, S. Ohnaka, J. da Providencia, and T. Kodama, 

Phys. Lett.B 226, 207 (1989). 

D.J. Thouless, Nucl. Phys. 22, 78 (1961). 

Y.K. Gambhir, P. Ring, and A. Thimet, Ann. Phys. (N.Y.) 198, 132 (1990). 

D.H. Youngblood, H.L. Clark, and Y.W. Lui, Phys. Rev. Lett. 82, 691 (1999). 

G.A. Lalazissis, J. K5nig, and P. Ring, Phys. Rev. C 55, 540 (1997). 

J.P. Blaizot, J.F. Berger, J. Decharge, and M. Girod, Nucl. Phys. A 591, 435 (1995). 

G. Colo, N. Van Giai, P.F. Bortignon, and M.R. Quagha, Phys. Lett. B 485, 362 (2000). 



19 



FIGURES 



o 

CO 



■116 



4 - 



> 

CD 



i2- 



CO 

O 



DC 







Sn 



RPA(full response) 



'll". 

I J . RPA(with Dirac states of a) 



RPA(without Dirac states) 

RPA(with Dirac states of co) 




NL3 



10 20 30 40 50 

E(MeV) 

FIG. 1. ISGMR strength distributions in ^^^Sn calculated with the NL3 effective interaction. 

The solid and long-dashed curves are the RRPA strengths with and without the inclusion of Dirac 
sea states, respectively. The dash-dot-dot (dash-dot) curve corresponds to calculations in which 
only vector mesons (scalar mesons) are included in the coupling between the Fermi sea and Dirac 
sea states. 
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FIG. 2. Effects of nuclear magnetism on the IS GMR strength distribution in ^^^Sn. Sohd curve: 
full RRPA calculation; dash-dotted curve: without the matrix elements of the spatial components 
of the vector meson fields; dashed curve: without the contribution of the Dirac sea to the matrix 
elements of the spatial components of the vector meson fields. The free Hartree response is also 
displayed. 
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